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ame of Document] SPECIFICATION 
[Title of the Invention] NOISE SUPPRESSING APPARATUS AND 

STORAGE MEDIUM 
[Scope of Demand for Patent] 

[Claim 1] A noise suppressing apparatus for subjecting an 
input image signal representing an input radiographic image to 
a noise suppression process for suppressing noise components 
included in the radiographic image, comprising: 

a band-limited-image-signal generation means for 
generating a plurality of band-limited image signals 
respectively representing a plurality of band-limited images 
belonging to a plurality of different frequency bands, based 
on the input image signal; 

index-value obtaining means for obtaining an index value 
indicating a degree of suppression of the noise component , 
based on information indicating an exposure dose with which the 
radiographic image has been produced; and 

noise suppressing means for processing each of the 
plurality of band-limited image signals so as to suppress noise 
in each of the plurality of band-limited images based on the 
obtained index value. 

[Claim 2] The noise suppressing apparatus as defined in 
Claim 1, 

wherein the index-value obtaining means obtains the 
index value for each of the band-limited image signals, and 

the noise suppressing means processes each of the 
plurality of band-limited image signals so as to suppress the 
noise in each of the plurality of band-limited images based on 
the index value obtained for each band-limited image signal. 

[Claim 3] The noise suppressing apparatus as defined in 
Claim 1 or 2, 

wherein the index-value obtaining means obtains the 
index value for each pixel of each of the plurality of 
band-limited images represented by the band-limited image 
signals, and 



the noise suppressing means processes each pixel so as 
to suppress the noise based on the index value obtained for each 
pixel, 

[Claim 4] The noise suppressing apparatus as defined in any 
5 one of claims 1 to 3, wherein the index-value obtaining means 
obtains a first evaluation value from a first one of the 
plurality of band-limited image signals belonging to a first 
one of the plurality of different frequency bands and a second 
evaluation value from a second one of the plurality of 

10 band-limited image signals belonging to a second one of the 
plurality of different frequency bands which is lower than the 
first one of the plurality of different frequency bands, 
determines weights based on the information indicating the 
exposure dose with which the radiographic image has been 

15 produced, for use in a weighted sum of the first and second 
evaluation values, obtains the weighted sum, and obtains based 
on the weighted sum the at least one index value indicating the 
degree of suppression of the noise for the first one of the 
plurality of band-limited image signals. 

20 [Claim 5] The noise suppressing apparatus as defined in 

Claim 4, wherein the index-value obtaining means obtains each 
of the evaluation values for each pixel of one of the plurality 
of band-limited images corresponding to each of the evaluation 
values, based on pixel values of the band-limited image in the 

25 vicinity of the pixel. 

[Claim 6] The noise suppressing apparatus as defined in 
Claim 4 or 5, wherein the index-value obtaining means obtains 
the index value based on a length and/or an orientation of the 
pixel vector using as an evaluation value a pixel vector at each 

30 pixel of the band-limited image represented by the band-limited 
image signal. 

[Claim 7] The noise suppressing apparatus as defined in 
Claim 6, wherein the index-value obtaining means obtains as the 
index value at least one of a degree of edge confidence, an 
35 amount of pixel energy, and a vector orientation. 
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[Claim 8] The noise suppressing apparatus as defined in any 
one of claims 1 to 7, wherein the noise suppressing means 
obtains a transformed image signal by subjecting the 
band-limited image signal to a predetermined process, obtains 
5 a weighted sum of the band-limited image signal and the 
transformed image signal by using weights determined based on 
the index value, and carries out the noise suppressing. 

[Claim 9] The noise suppressing apparatus as defined in 
Claim 8, wherein, the noise suppressing means obtains the 
10 transformed image signal by performing as the predetermined 
process the steps of: obtaining a pixel vector at each pixel 
of the band-limited image represented by the band-limited image 
signal; arranging an orientation-dependent filter based on a 
length and an orientation of the pixel vector; and making 
15 convolution of pixel values in the vicinity of each pixel in 
the band-limited images with the orientation-dependent filter. 

[Claim 10] A computer-readable storage medium storing a 
program which instructs a computer to execute on an input image 
signal representing an input radiographic image a noise 
20 suppression process for suppressing noise components included 
in the radiographic image, 

wherein the program instructs the computer to carry out 
the steps of: 

generating a plurality of band-limited image signals 
25 respectively representing a plurality of band-limited images 
belonging to a plurality of different frequency bands, based 
on the input image signal; 

obtaining an index value which indicates a degree of 
noise suppression based on the exposure dose with which the 
30 radiographic image has been produced; and 

processing each of the plurality of band-limited image 
signals so as to suppress noise in each of the plurality of 
band-limited images based on the obtained index value. 
[Detailed Description of the Invention] 
35 [0001] 
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[Technical Field of the Invention] 

The present invention relates to an apparatus for 
subjecting an input image signal representing a radio graphic 
image to a noise suppression process to suppress noise 
5 components included in the radiographic image, and a storage 
medium storing therein a processing sequence for carrying out 
the noise suppression process. 
[0002] 

[Description of the Related Art] 
10 Currently, when radiographic images obtained by using a 

computed radiography device (hereinafter referred to as a CR 
device) or the like are used in diagnosis, image processing such 
as frequency emphasis processing or gradation processing is 
performed on the obtained radiographic images before the 
15 radiographic images are displayed on CRT monitors as soft 
copies or recorded in films as hard copies, 
[0003] 

Radiographic images suffer from the problem that they 
tend to include noticeable quantization noise in low density 
20 areas corresponding to low-intensity radiation exposure. 
Therefore, various methods have been proposed for suppressing 
noise components included in image signals carrying 
radiographic images . 
[0004] 

25 For example, Japanese Unexamined Patent Publications No. 

6 ( 1994 ) -96200 proposes a method of obtaining an image with 
noise components therein being suppressed, wherein an image is 
transformed (decomposed) into a set of detail images 
(represented by band-limited image signals corresponding to 1 

30 to M resolution levels) ; a moving average of squared pixel 
values of each detail image in an N X N vicinity centered around 
each pixel of interest (i.e., a sum of the squared pixel values 
divided by N 2 ) is calculated as a local variance; a local 
variance corresponding to the peak in the histogram of the local 

35 variance is obtained as a noise variance; the local variance 



corresponding to each pixel is compared with the noise 
variance; when the local variance is comparable to or smaller 
than the noise variance, a portion of the band-limited image 
signal corresponding to the pixel is reduced; and thereafter, 
5 the set of detail images processed as above are inversely 
composed to the space of the original image by inverse 
multiresolution transformation so that an image in which the 
noise is suppressed is obtained. 
[0005] 

10 Japanese Unexamined Patent Publications No. 

6(1994)-96200 also discloses a technique for calculating noise 
variances for detail images at lower resolution levels based 
on a noise variance for a detail image at the highest resolution 
level (i.e., the finest-grained detail image). 

15 [0006] 

[Problems to be Solved by the Invention] 

According to the method disclosed in Japanese Unexamined 
Patent Publications No. 6 ( 1994 ) -274 615 mentioned above, noise 
is suppressed by using a variance obtained from variances 

20 (moving average) only of a relevant resolution signal (detail 
image) and a histogram of the local variances. Thus, at an area 
where noise level is low, such as an area directly exposed to 
X-rays, discrimination between the noise and edge is easy and 
the noise can be effectively suppressed. However, for the 

25 radiographed image exposed to a relatively small X-ray dose, 
for example, an area in which images of objects exist, not edge 
signals but noise signals increase and the histogram shape is 
different, as a result of which discrimination between the 
noise and edge becomes impossible. Consequently, when noise 

30 suppression is applied, edge information is suppressed 
together with the noise, which causes edge degradation and 
hence decreases sharpness of the image. 
[0007] 

In addition, even for images with identical noise levels, 
35 histograms obtained from radiographic images including images 
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of objects having complex structures and histograms obtained 
from radiographic images including images of objects having 
simple structures are different in shape. Therefore, it is 
difficult to favorably discriminate between edges and noise 
5 according to variations between objects. 
[0008] 

The present invention has been made in view of the points 
mentioned above. An object of the present invention is to 
provide an apparatus for effectively suppressing noise in or 

10 eliminating noise from a radiographic image regardless of an 
exposure dose with which the radiographic image has been 
produced, and reducing edge degradation which is likely to be 
caused by noise suppression in a noisy radiographic image. 
Another object of the present invention is to provide a 

15 computer-readable storage medium storing therein a processing 
procedure for conducting the noise suppression process. 
[0009] 

[Means for Solving the Problems] 

According to the present invention, there is provided a 

20 noise suppressing apparatus for subjecting an input image 
signal representing an input radiographic image to a noise 
suppression process for suppressing noise components included 
in the radiographic image, comprising: a 

band-limited- image-signal generation means for generating a 

25 plurality of band-limited image signals respectively 
representing a plurality of band-limited images belonging to 
a plurality of different frequency bands, based on the input 
image signal; index-value obtaining means for obtaining an 
index value indicating a degree of suppression of the noise 

30 component , based on information indicating an exposure dose 
with which the radiographic image has been produced; and noise 
suppressing means for processing each of the plurality of 
band-limited image signals so as to suppress noise in each of 
the plurality of band-limited images based on the obtained 

35 index value. 
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[0010] 

The information indicating an exposure dose with 
which the radiographic image has been produced may be 
information which directly indicates the exposure dose with 
5 which the radiographic image has been produced or may be 
information which indirectly indicates the exposure dose with 
which the radiographic image has been produced, for example, 
information indicating a menu of radiography, an age of a 
patient, a condition of radiography (e.g., a condition for 
10 irradiation by a radiographic apparatus), a normalization 
condition (see, for example, Japanese Unexamined Patent 
Publication No. 2 (1990) -108175) , or a pixel value (density 
value) of the radiographic image, 
[0011] 

15 The index value indicating a degree of suppression of the 

noise may be any value which can indicate a degree of 
suppression of the noise included in the radiographic image. 
For example, information indicating the exposure dose per se 
can be used as the index value. 

20 [0012] 

In the noise suppressing apparatus of the invention, it 
is preferable that the index-value obtaining means obtains the 
index value for each of the band-limited image signals, and the 
noise suppressing means processes each of the plurality of 

25 band-limited image signals so as to suppress the noise in each 
of the plurality of band-limited images based on the index value 
obtained for each band-limited image signal. 
[0013] 

Further, in the noise suppressing apparatus of the 
30 invention, it is preferable that the index-value obtaining 
means obtains the index value for each pixel of the band-limited 
images represented by the band-limited image signals, and the 
noise suppressing means processes each of the pixel so as to 
suppress the noise based on the index value obtained for each 
35 pixel. 
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[0014] 

Further, in the noise suppressing apparatus of the 
invention, it is preferable that the index-value obtaining 
means obtains a first evaluation value from a first one of the 
5 plurality of band-limited image signals belonging to a first 
one of the plurality of different frequency bands and a second 
evaluation value from a second one of the plurality of 
band-limited image signals belonging to a second one of the 
plurality of different frequency bands which is lower than the 

10 first one of the plurality of different frequency bands, 
determines weights based on the information indicating the 
exposure dose with which the radiographic image has been 
produced, for use in a weighted sum of the first and second 
evaluation values, obtains the weighted sum, and obtains based 

15 on the weighted sum the at least one index value indicating the 
degree of suppression of the noise for the first one of the 
plurality of band-limited image signals, 
[0015] 

In this case, it is preferable that the index-value 
20 obtaining means obtains each of the evaluation values for each 
pixel of one of the plurality of band-limited images 
corresponding to each of the evaluation values, based on pixel 
values of the band-limited image in the vicinity of the pixel. 
[0016] 

25 Alternatively, it is also preferable that the 

index-value obtaining means obtains the index value based on 
a length and/or an orientation of the pixel vector using as an 
evaluation value a pixel vector at each pixel of the 
band-limited image represented by the band-limited image 

30 signal. In this case, it is preferable that the index-value 
obtaining means obtains as the index value at least one of a 
degree of edge confidence, an amount of pixel energy, and a 
vector orientation . 
[0017] 

35 in the noise suppressing apparatus of the invention, it 
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is preferable that the noise suppression process is carried out 
by obtaining a transformed image signal by subjecting the 
band-limited image signal to a predetermined process, and 
obtaining a weighted sum of the band-limited image signal and 
5 the transformed image signal by using weights determined based 
on the index value. 
[0018] 

In this case, it is preferable that the noise suppressing 
means obtains the transformed image signal by performing as the 

10 predetermined process the steps of: obtaining a pixel vector 
at each pixel of the band-limited image represented by the 
band-limited image signal; arranging an orientation-dependent 
filter (anisotropic filter) based on a length and an 
orientation of the pixel vector; and making convolution of 

15 pixel values in the vicinity of each pixel in the band-limited 
images with the orientation-dependent filter. 
[0019] 

A storage medium of the invention instructs a computer 
to execute on an input image signal representing an input 

20 radiographic image a noise suppression process for suppressing 
noise components included in the radiographic image, wherein 
the program instructs the computer to carry out the steps of: 
generating a plurality of band-limited image signals 
respectively representing a plurality of band-limited images 

25 belonging to a plurality of different frequency bands, based 
on the input image signal; obtaining an index value which 
indicates a degree of noise suppression based on the exposure 
dose with which the radiographic image has been produced; and 
processing each of the plurality of band-limited image signals 

30 so as to suppress noise in each of the plurality of band-limited 
images based on the obtained index value. 
[0020] 

[Advantageous effect of the Invention] 

According to the noise suppression apparatus and the 
35 storage medium of the present invention, an index value 



indicating a degree of suppression of the noise is obtained 
based on information indicating an exposure dose with which the 
radiographic image has been produced, and each of a plurality 
of band-limited image signals is processed so as to suppress 
5 noise in each of the plurality of band-limited images based on 
the at least one index value. Therefore, when a radiographic 
image is reconstructed from the plurality of band-limited image 
signals after noise suppression or elimination is performed in 
each of the plurality of band-limited image signals, a 
10 radiographic image in which noise is effectively suppressed or 
eliminated can be obtained regardless of the exposure dose with 
which the original radiographic image has been produced, 
[0021] 

The index-value obtaining unit may obtain the at least 
15 one index value indicating the degree of suppression of the 
noise for each of the plurality of band-limited image signals, 
and the noise suppression unit may process each of the plurality 
of band-limited image signals so as to suppress the noise in 
each of the plurality of band-limited images based on the at 
20 least one index value obtained for one of the plurality of 
band-limited image signals representing the band-limited image. 
In this case, each frequency component of the input image signal 
can be processed by using at least one optimum index value 
obtained for each of the plurality of frequency bands. 
25 Therefore, fine-grained noise suppression can be performed. 
[0022] 

Further, when the index-value obtaining means obtains a 
first evaluation value from a first one of the plurality of 
band-limited image signals belonging to a first one of the 

30 plurality of different frequency bands and a second evaluation 
value from a second one of the plurality of band-limited image 
signals belonging to a second one of the plurality of different 
frequency bands which is lower than the first one of the 
plurality of different frequency bands, determines weights 

35 based on the information indicating the exposure dose with 
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which the radiographic image has been produced/ for use in a 
weighted sum of the first and second evaluation values, obtains 
the weighted sum, and obtains based on the weighted sum the at 
least one index value indicating the degree of suppression of 
5 the noise for the first one of the plurality of band-limited 
image signals, the processing for suppressing noise can be 
performed by using at least one proper index value which is 
determined in consideration of both of noise and edges. 
[0023] 

10 When the image is noisy, edge degradation is likely to 

occur even in an edge region by the noise suppression since the 
evaluation values are calculated based on the information of 
the band of interest involving a poor signal-to-noise ratio. 
However, when the first and second evaluation values are 

15 obtained based on pixel values for each pixel of the 
corresponding one of the plurality of band-limited images in 
the vicinity of the pixel, substantially, information in a 
lower frequency band in which the signal-to-noise ratio is 
increased is used for obtaining the first and second evaluation 

20 values. Therefore, the capability of discrimination between 
noise and edges is increased, and thus the aforementioned 
problem (edge degradation) can be reduced. 
[0024] 

Further, in this case, when the index value is obtained 
25 based on a length and/or an orientation of the pixel vector 
using as an evaluation value a pixel vector at each pixel of 
the band-limited image represented by the band-limited image 
signal, noise can be suppressed more effectively by using the 
at least one index value including information on 
30 characteristics (an orientation and sharpness) of the edge. 
Therefore, the edge degradation, which is likely to be caused 
by noise suppression in a noisy radiographic image, can be 
reduced. Further, variations in image quality due to variations 
in the exposure dose can be suppressed. 
35 [0025] 
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Further, when the index-value obtaining means obtains as 
the index value at least one of a degree of edge confidence, 
an amount of pixel energy, and a vector orientation, the noise 
on the edge can be properly suppressed since the band-limited 
5 image signals can be smoothed based on the vector information. 
Further, when all of the degree of edge confidence, the amount 
of pixel energy, and the orientation of the pixel vector are 
used, the noise on the edge can be further effectively 
suppressed without reducing edge contrast. 
10 [0026] 

Further, when the noise suppression process is carried 
out by obtaining a transformed image signal by subjecting the 
band-limited image signal to a predetermined process, and 
obtaining a weighted sum of the band-limited image signal and 
15 the transformed image signal by using weights determined based 
on the index value, it becomes possible to discriminate between 
edges and noise and to carry out the noise suppression processs, 
each corresponding to the noise and the edge. 
[0027] 

20 In this case, when the noise suppressing means obtains 

the transformed image signal by performing as the predetermined 
process the steps of: obtaining a pixel vector at each pixel 
of the band-limited image represented by the band-limited image 
signal; arranging an orientation-dependent filter 

25 (anisotropic filter) based on a length and an orientation of 
the pixel vector; and making convolution of pixel values in the 
vicinity of each pixel in the band-limited images with the 
orientation-dependent filter, it becomes possible to 
accurately discriminate between edges and noise, and the 

30 overall throughput of the processing can be increased. 
[0028] 

[Embodiments of the Invention] 

Hereinafter, embodiments of the present invention will 
be described in detail with reference to the drawings. 
35 [0029] 
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Fig. 1 is a diagram illustrating an outline of a 
construction of a noise suppressing apparatus according to the 
present invention. 

[0030] 

5 As shown in FIG. 1, a noise suppressing apparatus 100 

comprises: a band-limited- image-signal generation unit 1, the 
band-limited- image-signal generation unit 1 generating a 
plurality of band-limited image signals that represent a 
plurality of band-limited images respectively belong to a 

10 plurality of different frequency bands, based on an input image 
signal Sin that represents a radiographic image obtained by an 
image reading device or the like and having a certain 
resolution; an index-value obtaining unit 2, the index-value 
obtaining unit 2 obtaining at least one index value which 

15 indicates a degree of noise suppression based on information 
indicating an exposure dose with which the radiographic image 
has been produced; a noise-suppression processing unit 3, the 
noise-suppression processing unit 3 performing noise 
suppression processes on each of the plurality of band-limited 

20 image signals according to the degrees of noise suppression; 
and an image reconstruction unit 4, the image reconstruction 
unit 4 composing (reconstructing) a processed image signal 
Sproc which represents a noise-suppressed radiographic image, 
from the plurality of band-limited image signals on which the 

25 noise suppression process is performed by the noise 
suppression processing unit 3. 
[0031] 

Further, according to the particular embodiment, for 
example, in a radiographic-image-inf ormation 

30 recording-and-reproducing system as disclosed in Japanese 
Unexamined Patent Publication Nos . 55 ( 1980) -12492 and 
56 ( 1981) -11395, a radiographic image of a human body that is 
recorded in a stimulable phosphor sheet and read as a digital 
image signal by laser beam scanning is targeted for the noise 

35 suppression process. The radiographic image is read by 
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performing the laser beam scanning in two dimensions by moving 
a laser beam on the stimulable phosphor sheet in a main scanning 
(lateral) direction while moving the stimulable phosphor sheet 
in a feeding (longitudinal) direction. 
5. [0032] 

Operations of a noise suppressing apparatus 100 having 
the structure described above are explained below. 
[0033] 

First, an outline of the processing is described with 
10 reference to a flowchart shown in FIG. 2. 
[0034] 

In order to generate the plurality of band-limited image 
signals, it is preferable to use the multiresolution 
transformation such as the Laplacian pyramid decomposition 

15 that is proposed in Japanese Unexamined Patent Publication Nos. 
5 (1993) -244508 and 6 ( 1994 ) -0962 00 and Japanese Patent 
Application Nos. 11 (1999) -363766 and 2000-022828 or the 
wavelet transformation that is proposed in Japanese Unexamined 
Patent Publication No. 6 ( 1994 ) -274 615 and Japanese Patent 

20 Application No. 11-363766 which are also assigned to the 
assignee of the present invention. Alternatively, the 
plurality of band-limited image signals may be generated by 
using other known methods. For example, the plurality of 
band-limited image signals may be generated by using the 

25 unsharp mask signals as disclosed Japanese Unexamined Patent 
Publication No. 10 (1998) -75364 . Descriptions on embodiments 
below are based on the assumption that the Laplacian pyramid 
decomposition is used. 
[0035] 

30 Band-limited image signals are generated from an entered 

original image by the use of the Laplacian pyramid 
decomposition that is one of the multiresolution 
transformations (step S21). Then, vector components at each 
pixel position of each of a plurality of band-limited images 

35 in multiresolution spaces are calculated, where the 
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band-limited images are respectively represented by the 
band-limited image signals (step S22) . When the vector 
components are obtained in double-angle representation (which 
is explained later) , four-orientation vector components, i.e., 
5 vector components corresponding to four orientations at 
intervals of 45 degrees, are obtained in each pixel position. 
Based on the four-orientation vector components, it is possible 
to discriminate between noise components and edge components 
in each pixel position. 

10 [0036] 

When a singular point (local noise) exists, e.g., when 
a vector which is extremely greater than surrounding vectors 
exists, the local noise at each pixel position is likely to be 
incorrectly recognized as an edge signal. Therefore, the 

15 vicinity average (vector average), i.e., an average of values 
of each vector component in the vicinity, is obtained by using 
a one-dimensional filter (step S23) . This operation is based 
on an assumption that edge signals are continuous. In this 
embodiment, the vicinity average is obtained by using an 

20 isotropic two-dimensional space filter. Further, the vector 
average is modified in by using a vector component obtained 
at a resolution level lower than the resolution level at which 
the vicinity average is obtained. At this time, the 
modification of the vector average is made based on information 

25 on the exposure dose with which the original radiographic image 
has been produced (step S24) . 
[0037] 

Next, a degree C of edge confidence and an index E of pixel 
energy are calculated based on each vector averaged and 

30 modified as above in accordance with a method explained later 
(step S25) . Noise suppression process using adaptive filtering 
is performed based on the degree C of edge confidence and the 
index E of pixel energy (step S26) . Finally, Laplacian pyramid 
reconstruction that is one of the inverse multiresolution 

35 transformations is made so that a processed image in which noise 
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is suppressed is obtained (step S27). 
[0038] 

The adaptive filtering in step S26 is .performed by an 
anisotropic filter (orientation-dependent filter) and an 
5 isotropic filter (orientation-independent filter) . It is 
possible, for example, to calculate tens of different sets of 
anisotropic filter coefficients for the anisotropic filter in 
advance, and select one of the tens of different sets according 
to a vector orientation D. On the other hand, the isotropic 
10 filter can be realized by a simple non-linear transformation. 
[0039] 

Next, details of the operations performed in the sequence 
are explained below. 
[0040] 

15 Fig. 3 is a block diagram illustrating an outline of the 

construction of the band-limited- image-signal generation unit 
1 and schematically illustrating the operations of generating 
five band-limited image signals corresponding to five 
resolution levels . 

20 [0041] 

For example, as disclosed in Japanese Unexamined Patent 
Publication No . 5 ( 1993 ) -244508 , the filtering processing unit 
10 performs filtering processing on the input image signal Sin 
in each of the main scanning and feeding directions so as to 

25 produce a signal LI (hereinafter referred to as a lower 
resolution signal) , which belongs to a lower resolution level 
than the resolution level of the input image signal Sin. Then, 
the filtering processing unit 10 performs filtering processing 
on the lower resolution signal LI in each of the main scanning 

30 and feeding directions so as to produce a next lower resolution 
signal L2. Thereafter, further lower resolution signals Lk (k=l 
to n) are successively obtained by repeating the filtering 
processing in a similar manner. Next, the interpolation 
processing unit 11 performs interpolation processing on each 

35 of the lower resolution signals Lk so that the number of pixels 



16 



in each of the main scanning and feeding directions is doubled, 
i.e., the number of pixels in the lower resolution signal Lk 
is increased by a factor of four. Thus, a plurality of unsharp 
image signals Susl to Susn (hereinafter collectively referred 
5 to as "Susk (k=l to n)") each having a different degree of 
sharpness are obtained. Thereafter, the subtractor 12 obtains 
a difference between each of the lower resolution signal Lk-1 
and one of the plurality of unsharp image signals Susk having 
the same number of pixels as the lower resolution signal Lk-1 
10 so as to generate one of the band-limited image signals Bk. 
[0042] 

Next, details of the operations performed for obtaining 
an index value (a degree of noise suppression) and suppressing 
noise by using the band-limited image signals Bk obtained as 
15 described above are explained below. 
[0043] 

Fig. 5 is a block diagram illustrating details of a 
general configuration of the noise suppressing apparatus 100. 
As shown, the index-value obtaining unit 2 comprises for each 

20 of the plurality of band-limited image signals Bk: a 
pixel-vector generation unit 22 for generating a pixel vector 
at each pixel of a band-limited image represented by the 
corresponding band-limited image signal Bk; and an index-value 
calculation unit 24 for obtaining at least one of a degree of 

25 edge confidence, an index of pixel energy, and an edge 
orientation (which are examples of the index values of the 
present invention) , for each pixel of the band-limited image 
represented by the corresponding band-limited image signal Bk, 
based on the length and/or orientation of the pixel vector 

30 generated by the pixel-vector generation unit 22 corresponding 
to the same resolution level. 
[0044] 

The noise-suppression processing unit 3 comprises a 
suppression processing unit 32 for each of the plurality of 
35 band-limited image signals Bk. The suppression processing unit 
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32 serves to suppress noise included in the band-limited image 
signal Bk, based on the index value output from the index-value 
calculation unit 24. 
[0045] 

5 The noise suppression process of the invention is based 

on a technique that "a line signal is smoothed along the 
orientation of the line, and isolated noise is 
two-dimensionally smoothed" . The most characteristic feature 
of the noise suppression process is that a smooth edge is 

10 obtained by the smoothing operation of the line (edge) signal, 
and information necessary for the smoothing operation is 
represented in only a vector or tensor form. In the examples 
disclosed in this specification, as shown in FIG, 2, the 
double-angle representation (hereinafter referred to as "D-A 

15 representation") is used as a vector representation form. 
[0046] 

The D-A representation of a vector is a technique for 
representing a line signal, and advantageous in that the degree 
of confidence of a line signal (index of lineness) can be 
20 obtained by only calculating the vicinity average of 
information in the D-A representation. This feature is 
explained below with reference to FIG. 6. 
[0047] 

When density vectors of an image signal as shown in FIG. 

25 6(A) are calculated, the density vectors are represented as 
illustrated in FIG. 6(B) in the full-angle representation 
(hereinafter referred to as "F-A representation") , which is an 
example of normal vector representations. Thus, the 
orientations of the vectors on the respective sides of the 

30 boundary are opposite, where the low-density area in the image 
corresponds to the boundary. On the other hand, when 
represented in the D-A representation, the orientations of the 
vectors on both sides of the boundary are identical as 
illustrated in FIG. 6(C) since the angle values of the vectors 

35 are doubled in the D-A representation. 
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[0048] 

In addition, when the degrees C of edge confidence of the 
above. vectors are obtained by the vicinity average as indicated 
with the bold arrows in FIGS. 6(B) and 6(C), the degree C of 
5 edge confidence shown in FIG. 6 (C) is considerably smaller than 
the degree C of edge confidence shown in FIG. 6(B) . Although 
not specifically shown, it will be easily understood that the 
confidence of noise also becomes small (since vectors 
surrounding a vector of interest have random orientations) . 
10 Therefore, it is difficult to discriminate between noise and 
line information in the F-A representation. 
[0049] 

On the other hand, in the D-A representation, vectors 
indicating line (edge) orientations can be defined as 

15 illustrated in FIG. 7. In the drawing, qO to q3 each indicate 
the magnitude of a directional component at a pixel of interest. 
When the magnitudes of two orthogonal directional components 
at a pixel are identical (i.e., when the pixel is located at 
a point of intersection of two lines) , the output in the D-A 

20 representation becomes small. When two orthogonal directional 
components have different magnitudes, the orientation 
corresponding to the greater magnitude becomes a main 
orientation at the pixel. 
[0050] 

25 Therefore, when the four directional components qO to q3 

of a vector at each pixel are obtained, the vector can be 
represented in the D-A representation. 
[0051] 

Details of operations for obtaining the directional 
30 components qO to q3 are described below. 
[0052] 

The band-limited image signals used for calculation are 
Laplacian signals generated by the Laplacian pyramid 
decomposition. The four directional components are calculated 
35 by convolution of the Laplacian signals with the four 
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two-dimensional filters as illustrated in FIG. 8. Table 4-1 
shows examples of filter coefficients of a 5 X 5 qO filter. Since 
each of .the Laplacian signals can become positive or negative, 
and each of the filter coefficients can be positive or negative, 
5 absolute values of the convolution products are used for 
calculation of the directional components. 
[TABLE 1] 



Filter Coefficients for qO Filter (5-by-5) 



0.0012 


0.0211 


0. 0577 


0 . 0211 


0.0012 


0.0053 


0 .1389 


0. 6093 


0 . 1389 


0.0053 


0.0000 


0.0000 


0.0000 


-0. 0000 


-0.0000 


-0.0053 


-0.1389 


-0. 6093 


-0. 1389 


-0.0053 


-0.0012 


-0.0211 


-0. 0577 


-0.0211 


-0.0012 



10 [0053] 

Each of the above-mentioned four filters is a kind of 
differential filter. Therefore, in the case where a 
second-derivative signal such as a Laplacian signal is 
convoluted, the output of the filter does not become great in 

15 regions in which the gradient of the second-derivative signal 
is not great, even when the second-derivative signal per se is 
great in the regions. This feature is explained below with 
reference to FIG. 9. 
[0054] 

20 FIG. 9 is a diagram illustrating a relationship between 

a Laplacian signal and an output of a first-derivative filter 
(the absolute value of the first derivative of the Laplacian 
signal) . As shown, the point a of the Laplacian signal 
corresponds to a portion of an edge region, and the Laplacian 

25 signal is maximized at the point A. However, since the gradient 
of the Laplacian signal at the point A is zero, the output of 
the first derivative corresponding to the point a becomes zero. 
[0055] 

In addition, at a point C which is a boundary point 
30 between the edge region and a non-edge region, the output of 
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the first derivative corresponding to the point C becomes 
grater than the output of the first derivative corresponding 
to the point A. This tendency is strengthened with increase of 
the mask size. However, when the mask size is increased, 
5 trackability of very small edges decreases, and therefore 
sharpness of images of the very small edges decreases, 
[0056] 

Thus, a small mask size is preferable when the image 
signal does not include artificially produced noise as in a 
10 pattern image. In an actual input image as it is, however, noise 
exists in the entire image in varying degrees. Therefore, first 
derivatives generated by using filters having a small mask size 
are strongly affected by noise. 
[0057] 

15 The above circumstances can be a cause of image quality 

degradation when an adaptive filter which operates based on 
pixel energy (i.e., an average of the directional components 
qO to q3) is used. 
[0058] 

20 Nevertheless, when the amount of noise can be estimated, 

it is possible to optimize the filter by changing the mask size 
according to the amount of noise, or setting filter 
coefficients so as to substantially achieve the effect of the 
setting of the mask size. For example, when a radiographic image 

25 is to be processed as in the present embodiment, an X-ray dose 
and an amount of noise can be estimated from the S value 
(indicating the reading sensitivity) and the L value 
(indicating the latitude) , and an optimum set of filter 
coefficients can be calculated. Regarding the S value and the 

30 L value, see, for example, Japanese Unexamined Patent 
Publication No. 2 (1990) -108175 . 
[0059] 

Specifically, the vicinity average of each of the vector 
components (qO to q3) at each pixel is obtained. For the 
35 vicinity average, an isotropic two-dimensional filters as 
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illustrated in FIG. 10 is used. 
[0060] 

When the mask size of the two-dimensional filter is 
varied, the smoothing level of the vector component of course 
5 varies. The smoothing level is reflected on the edge confidence 
and the pixel energy, and the influence of the smoothing level 
on the final image is relatively great. When the mask size is 
increased, noise and relatively large edges can be 
discriminated with high accuracy, and small edges are likely 

10 to be regarded as noise. Therefore, smoothing with a large mask 
size is effective when the input image does not include fine 
structures, for example, as in a radiographic image of a chest 
of a child. On the other hand, since bone images such as an image 
of a foot includes complex fine structures such as trabecula 

15 here and there, fine signals indicating the complex fine 
structures cannot be recognized when the smoothing level is 
raised. Therefore, a small mask size is used. 
[0061] 

As proposed in Japanese Patent Application No. 

20 2000-022828 by the present applicant, the above vicinity 
average can be modified by using vector components in a lower 
resolution image which has a lower resolution than the targeted 
image. Investigations by the present inventors have revealed 
the following characteristics when the vector calculated from 

25 a currently targeted band-limited image signal (band-limited 
image of interest) for vector averaging and another vector 
(having a lower resolution) are used. 

1) An image including therein a smaller amount of noise 
(the image which undergoes a large amount of X-ray exposure) 

30 has a higher signal-to-noise ratio. Therefore, vector 
components obtained by the use of the vector average of the 
currently targeted band-limited image signals can more 
faithfully follow fine signals representing fine structures in 
the input image than vector components modified with the lower 

35 resolution signals, which prevents the edge degradation. 
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2) An image including therein a larger amount of noise 
has a lower signal-to-noise ration. Therefore, when an average 
of a vector .calculated from a band-limited image signal in a 
frequency band (with a lower SNR) of interest and another vector 
5 calculated based on image information in a lower-frequency band 
(with an enhanced SNR) is used, the average of the vectors can 
more faithfully follow relatively large signals which are not 
buried in the noise in the input image, and noise can be 
effectively suppressed. 
10 Thus, in order to improve image quality, it is important to 
control the degree of noise suppression according to the amount 
of noise included in the input image. 
[0062] 

Then, a noise amount evaluating process which becomes 
15 necessary for controlling the degree of noise suppression is 
described below. 
[0063] 

For a radiographic image such as an X-ray photogram, 
noises are mainly caused by reduction of transmission radiation 
20 dose. Therefore, if the transmission radiation dose is known, 
an approximate noise amount could be estimated. 
[0064] 

The average C of vectors having respectively different 
frequency bands can be calculated in accordance with the 
25 equation (1) given below, where A is a average of vectors 
(vector average) calculated from a band-limited image signal 
in a frequency band of interest, B is another average of vectors 
calculated based on image information in a lower frequency 
band. 
3 0 [Formula 1] 

C = f(x)xA + (l-f(x))xB (1) 
[0065] 

In the equation, x represents an X-ray dose, and f (x) is 
a function of the X-ray dose x, and represents a weight of the 
35 vicinity-averaged vector in the weighted average. 
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[0066] 

For estimating the amount of noise, various information 
which may indicate X-ray dose can be used, for example,. 1) the 
radiographed region or the menu of radiography; 2) the S or L 
5 value indicating a normalization condition (EDR condition) ; 3) 
the signal values (density values) of the image; and 4) the age 
of a patient or the condition of radiography, 
[0067] 

When the radiographed region or the menu of radiography 
10 is used, for example, a low-dose menu, a child menu, and the 
like may be provided, and the weights of the aforementioned 
vicinity- averaged vectors in the weighted average can be 
changed according to a selected one of the menu. When the S or 
L value indicating a normalization condition is used, it is 
15 possible to use the function f (x) which increases with decrease 
in the S value (corresponding to increase in the X-ray dose) . 
For example, it is preferable to use, for example, the functions 
f (x) given by the following set of equations (2) . It is of course 
that other such functions may be used. 
2 0 [Formula 2] 

f(x) = 1.0 when S < 100 

f(x) = (2000 - S) / 1900 when 100 < S < 200ol (2) 

f(x) = 0.0 when S > 2000 \ 

[0068] 

When the signal value of the image is used, since the 
density value corresponds to the X-ray exposure on a stimulable 
25 phosphor sheet, the relative amount of the X-ray exposure dose 
can be indicated by using the signal value. Therefore, it is 
preferable to use the functions, for example, as given by the 
following set of equations (3) . 
[Formula 3] 

x = s x io ( - LxQL/ 1024) ' 
f(x) = 1.0 when x < 100 

f(x) = (2000 - x) / 1900 when 100 < x < 2000 
f(x) = 0.0 when x > 2 000 



30 



(3) 
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wherein QL denotes a signal value, and x denotes a relative 
amount of the X-ray dose. 
[0069] . 

Alternatively, the set of equations (3) may be modified 
5 so that the function f (x) satisfies 0 . 5<f (x) <1 . 0, instead of 
0<f (x)<1.0. 

[0070] 

That is, in the case where the signal values are used, 
the signal value of each pixel is referred to when the vector 

10 average is calculated, the relative amount of the X-ray 
exposure dose is estimated in accordance with the set of 
equations (3) , and a weighted average of vectors is calculated 
by using the weights of the vicinity-averaged vectors 
determined in accordance with the function defined in advance. 

15 [0071] 

Next, the orientations and lengths of a primary vector 
and a secondary vector are calculated by using the four vector 
components (qO to q3) . Since each vector component is obtained 
in the D-A representation as illustrated in FIG. 7, the length 

20 VI of the primary vector is calculated in accordance with 
equations (4), and the unit-vector components exl and eyl of 
the primary vector are calculated in accordance with the 
equations (5) . 
[Formula 4] 

Zl = q0-q2 

25 Z2 = ql-q3 I (4) 

Vl = ( z l 2 +z2 2 ) 1/2 j 
[Formula 5] 

exl = zl / VI] 
eyl = z2 / VI J 
[0072] 

Since the orientation of the secondary vector is opposite 
30 to that of the primary vector in the D-A representation, the 
unit-vector components (ex2, ey2) of the secondary vector 
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oriented orthogonal to the primary vector are calculated in 
accordance with equations (6). 
[Formula 6] 

ex2 = -exl ] 

y (6) 

ey2 = -eyl J 
5 [0073] 

In addition, the amount of pixel energy Ve at each pixel 
is defined as an average of the components in accordance with 
the equations (7) . Further, the length V2 of the secondary 
vector can be calculated from the ratio of the pixel energy Ve 
10 and the length VI of the primary vector in accordance with the 
equations (8) 
[Formula 7] 

Ve= (q0+ql+q2+q3) /4 (7) 
[Formula 8] 

15 V2= (1-Vl/Ve) xVl (8) 

[0074] 

Next, by using the two pieces of information VI and V2 
on the mutually orthogonal primary and secondary vectors which 
have been averaged and modified based on the noise amount, the 

20 degree C of edge confidence, the index E of pixel energy, and 
a smoothing direction D of the anisotropic filter are 
calculated in accordance with the following equations. That is, 
the index E of pixel energy is calculated in accordance with 
the set of equations (9) by using a predetermined threshold 

25 value Th, and the degree C of edge confidence is calculated from 
the lengths VI and V2 of the primary and secondary vectors in 
accordance with the equation (10). 
[Formula 9] 



E = (Ve / Th) 2 / 2 when Ve < Th 

E = {1 - (2 - Ve / Th) 2 } / 2 when Th < Ve < 2 x Th 
E = 1.0 when 2 x Th < Ve 

30 [Formula 10] 

C= (V1-V2) /VI (10) 



(9) 
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[0076] 

Note that, 0.0 < E < 1.0, 0.0 £ C £ 1.0, 0.0 <D< 
31 (D is an integer) , 
[0075] 

5 Note that, 0.0<E<1.0, 0.0<C<1.0, 0.0<D<31 (D is an 

integer) . 

[0076] 

Further, as shown in the set of equations (11) , an angle 
9 is calculated from the unit-vector components of the 
10 secondary vector, and a quantized smoothing direction D of the 
anisotropic filter is calculated from the calculated angle. In 
this particular embodiment, 32 discrete angle values are 
employed. 

[Formula 11] 



6 - cos 1 (ex2) when ey2 > 0 
15 0 = cos -1 (-e*2) wheney2<0 

D = f(0) 



(13) 



wherin f (0 ) is the function that converts the continuous 
angle into discrete angle values. 
[0077] 

The degree C of edge confidence indicates a degree of 
20 likelihood that the pixel constitutes an edge, and becomes 
larger as the relevant pixel is positioned closer to a line and 
the larger index E of pixel energy indicates a degree of 
likelihood that the pixel constitutes a signal. A recognition 
model (an example of an adaptive filter) of a line, a point of 
25 intersection, an end point, and noise based on the degree C of 
edge confidence and the index E of pixel energy are indicated 
in FIG. 11. 

[0078] 

Next, smoothing processing using an anisotropic filter 
30 and adaptive filtering are performed on each band-limited image 
signal (Laplacian signal) at each pixel based on the degree C 
of edge confidence, the index E of pixel energy, and the 
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smoothing direction D, which are obtained as above from the 
information on the two mutually orthogonal vectors. 
[0079] 

In the present embodiment, smoothing processing is 
5 performed by using an anisotropic filter which is oriented 
along the line recognized based on the degree C of edge 
confidence since the larger the degree C of edge confidence, 
the greater the likelihood that a line is constituted. That is, 
the two-dimensional anisotropic spatial filter smoothes the 

10 Laplacian signal along the orientation of the primary vector. 
This anisotropic filter involves coefficients as illustrated 
in FIG. 12. Since noise in a noisy image is also superimposed 
on edge signals, the anisotropic filter as above is used for 
suppressing noise on an edge without reducing edge contrast. 

15 [0080] 

Next, an anisotropic filter as an orientation-dependent 
filter is selected depending on the smoothing direction D(0) 
which is calculated based on vectors. Then, the band-limited 
image signal is convoluted with the selected filter so as to 

20 produce a convolution product as an anisotropic-f iltered 
signal. The convolution product (anisotropic-f iltered signal) 
is designated by A in equation 12 and equation 13. 
[0081] 

FIG. 12 shows 5 segments an anisotropic filter in the 
25 range of 0 to 180 degrees which are obtained by dividing the 
360 degree filter into 8 segments. Further, although 
anisotropic filters preserve edges, the anisotropic filters 
smooth points at which edges intersect. However, intersection 
points appearing on actual images are not ideal intersection 
30 points (wherein sharp lines intersect at the right angle) , and 
therefore not so smoothed. 
[0082] 

As adaptive filtering for suppressing the noise 
components, the noise-suppression processing unit 3 calculates, 
35 as described above, a convolution product A of the 
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anisotropic-f iltered signal and the Laplacian signal and 
controls the weight of the convolution product A and the 
Laplacian signal based on the index E of pixel energy and the 
degree C of edge confidence so as to obtain a processed 
5 band-limited image signal fBk (k=l to n) (represented by Proc 
in the equations) for each pixel of the band-limited image 
represented by the band-limited image signal (Laplacian 
signal) . The processed band-limited image signal Proc is a 
signal in which noise components are suppressed. Definitions 
10 applied when calculating the band-limited image signal Proc are 
Definitions 1 and 2 described below. 
[0083] 

Definition 1: When FIG. 11(b) is defined as a linear: 
This is the case where FIG. 11(b) is defined as a linear, and 
15 the processed band-limited image signal Proc in which noise 
components are suppressed is calculated in accordance with the 
following equation (12) . 



When the degree C of edge confidence is high, an 
anisotropic filter output is selected because of edges (first 
term) , while when the degree C of edge confidence is low, the 
original Laplacian signal is attenuated by the index E of pixel 
25 energy. Thus, noise and intersection points are separated from 
edges (second term) . 



Definition 2: When FIG. 11(b) is defined as a noise: 
This is the case where FIG. 11(b) is defined as a mpose, and 
30 the processed band-limited image signal Proc in which noise 
components are suppressed is calculated in accordance with the 
following equation (13). 



20 



[Formula 12] 

Proc=CxA + Ex(l-C)xOrg , 
[0084] 



(12) 



[0085] 



35 



[Formula 13] 

Proc=£xCx^+£x(l-C)xO/-g . 
[0086] 



(13) 
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When the degree C of edge confidence is high, the 
anisotropic-f iltered signal is attenuated by the index E of 
pixel energy, so that the noise and edges are separated from 
each other (first term) . When the degree C of edge confidence 
5 is low, the Laplacian signal is attenuated by the index E of 
pixel energy, so that noise and intersection points are 
separated from each other (second term) . 
[0087] 

However, in the case of Definition 1, noise can be 
10 eliminated from the processed band-limited image signal fBk 
(Proc) only when the signal-to-noise ratio (SNR) is high. 
However, when the signal-to-noise ratio is low, noise 
elimination effect is not substantially produced. From these 
results, it is considered the degree of edge confidence at a 
15 pixel can become high in the low-SNR situation even when the 
pixel represents noise. 
[0088] 

In the case of definition 2, an artifact resulting from 
discontinuity (local discontinuity caused by noise or an 

20 artificially produced edge which suddenly appears) is more 
likely to be produced as the signal-to-noise ratio decreases. 
In order to eliminate the artifact, the threshold value used 
in the calculation of pixel energy is required to be raised. 
However, when the threshold value is raised, edges become 

25 unclear. This is because the pixel energy of edge competes with 
the pixel energy of noise in the low-SNR situation. This problem 
is inevitable as long as the pixel energy is used. 
[0089] 

Further, the degree C of edge confidence of a pixel has 
30 a small value when the pixel constitutes a nonlinear signal 
(representing an intersection point, an end point, or the like) . 
Therefore, in order to discriminate between the nonlinear 
signal and noise, it is possible to determine whether each pixel 
corresponds to a true signal or noise by comparing the index 
35 E of pixel energy with a predetermined value. In practice, it 
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is possible to continuously determine the unlikelihood of noise 
by using an arbitrary nonlinear function. For example, the 
unlikelihood N of noise can be determined by using a nonlinear 
function of a threshold value TH and the index E of pixel energy, 
5 as indicated in the equation (14). 
[Formula 14] 

exp(TH/E)-l x2x JL (14) 
exp (TH / E) + 1 TH 

[0090] 

In this case, it is preferable to replace the 
10 aforementioned equation (12) with the equation (15), and the 
aforementioned equation (13) with the equation (16) based on 
the degree C of edge confidence and the unlikelihood N of noise . 
[Formula 15] 

Proc=C xA + Nx(l-C)xOrg (15) 
15 [Formula 16] 

Proc=ExCxA + Nx(l-C)xOrg (16) 
[0091] 

Thus, after the noise suppression process is performed 
on the band-limited image signals, i.e., the processed 

20 band-limited image signals Proc are obtained, the image 
reconstruction unit 4 makes a Laplacian reconstruction as an 
inverse multiresolution transformation so as to produce a 
processed image signal Sproc which represents an image in which 
noise components are suppressed. 

25 [0092] 

As illustrated in FIG. 5, the image reconstruction unit 
4 comprises an interpolation processing unit 43 and an adder 
44 corresponding to each resolution level. The interpolation 
processing unit 43 performs interpolation processing on the 

30 processed band-limited image signal. The adder 44 obtains a sum 
of the processed band-limited image signal and the interpolated 
(magnified) image signal. 
[0093] 

FIG. 13 is a diagram schematically illustrating an 
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operation for performing the Laplacian reconstruction. After 
the processed band-limited image signals fBk in which noise 
components are suppressed . (k=l to n) are obtained, the 
interpolation processing unit 43 performs interpolation 
5 processing, as with the aforementioned interpolation 
processing unit 11, on the processed band-limited image signal 
fBn at the lowest resolution level so as to obtain an 
interpolated magnified image signal Sn f . Then, the magnified 
noise component signal Sn' is added to the noise component 

10 signal fBk-1 by the adder 44 so as to produce an added noise 
component signal Sn-1 at the second lowest resolution level. 
Such processing is repeatedly performed so as to obtain a higher 
resolution, thereby providing an added noise component signal 
SI with the highest resolution level. The added noise component 

15 signal SI with the highest resolution level is processed so as 
to obtain a noise-suppressed image signal Sproc. 
[0094] 

Thus, when an image is output based on the processed image 
signal Sproc, "an image with a low exposure dose can obtained 

20 an image having an appearance which is similar to the 
appearances of slightly unsharped images of high-frequency 
components" by controlling parameters in the adaptive filter 
and the like. Since the weights in the weighted average of the 
vicinity averages of vectors at a resolution level of interest 

25 and a lower resolution level are controlled based on the X-ray 
dose or doses in the input image, fine-edge-oriented noise 
suppression is performed in areas of an image which are exposed 
with a high X-ray dose, and large-edge-oriented noise 
suppression is performed in areas of the image which are exposed 

30 with a low X-ray dose. Therefore, noise can be effectively 
suppressed or eliminated in the entire image, edge degradation, 
which can be caused by noise suppression in a noisy image, can 
be reduced, and variations in image quality caused by 
variations in the exposure dose can be suppressed. In other 

35 words, even when the amount of noise included in the image 
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varies due to variations in the exposure dose, the noise 
suppressing apparatus can effectively suppress the noise in the 
image, reduce artifacts (unnaturalness ) such as an arabesque 
pattern, and make the image more natural. In addition, 
5 degradation of fine signals can be reduced. Thus, it is possible 
to obtain high quality images. 
[0095] 

While the noise suppressing apparatus of the invention 
has been described in terms of preferred embodiments, it should 
10 be understood that the present invention not necessarily 
limited thereto. 
[0096] 

Although in the described embodiments, the band-limited 
image signals in respectively different frequency bands are 
15 generated from the input image signal Sin by the Laplacian 
pyramid decomposition, the band-limited image signals may be 
generated by the wavelet transformation as disclosed in 
Japanese Unexamined Patent Publication No. 6 ( 1 994 ) -274 615 . 
[0097] 

20 Although in the described embodiments, the 

noise-suppressed image signal is obtained from the plurality 
of noise-suppressed band-limited image signals by the inverse 
multiresolution transformation, an added noise component 
signal SHI representing all of the noise components may be 

25 obtained by using the multiple resolution signals, and 
subtracted from the input image signal Sin so as to obtain the 
noise-suppressed image signal Sin, as disclosed in Japanese 
Patent Application No. 11 ( 1999) -3637 66 . Fig. 14 is a block 
diagram illustrating details of such a noise suppressing 

30 apparatus 100. 

[0098] 

In the described embodiments, vector information is used 
as information locally calculated from pixel values in the 
vicinity of a pixel of interest (hereinafter referred to as 
35 local information) , and the index values indicating the degree 
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of noise suppression is obtained based on the vector 
information (the weighted average of vicinity averages of the 
vectors) , and filter processing is performed based on the index 
values. Alternatively, a tensor average or a moving average 
5 deviation may be used as the local information, instead of the 
above vector information. In the case where the moving average 
deviation is used, when a first moving average deviation 
corresponding to a first band-limited image signal is 
calculated, a second moving average deviation may be added to 

10 the first moving average deviation according to information 
which indicates an exposure dose such as a menu or condition 
of radiography or an amount corresponding to the radiation dose, 
where the second moving average deviation is calculated for the 
pixel of interest from a second band-limited image signal at 

15 a second resolution level which is lower than the first 
resolution level. 
[0099] 

In the described embodiments, the index values 
indicating the degree of noise suppression are obtained from 

20 the image signal (particularly from the band-limited image 
signals) , and the noise suppression process is performed based 
on the index values so as to obtain an image in which noise is 
effectively suppressed regardless of the exposure dose. 
However, the index values which indicate the degree of 

25 suppression of noise components may be defined in any other way. 
For example, information indicating the exposure dose with 
which the input image has been produced can be used as the index 
values. In this case, the information indicating the exposure 
dose with which the input image has been produced may be 

30 obtained, for example, by using a photo timer or the like. 
[0100] 

Further, for example, as disclosed in Japanese 
Unexamined Patent Publication No. 6 ( 1994 ) -96200, when the 
local variance of a detail image within an area that is lower 
35 than that when the local variances for respective detail images 
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are calculated, a weight of the second local variance in the 
addition may be changed according to the exposure dose at each 
location of the image, and an unlikelihood N of noise may be 
calculated so as to transform the original signal (band-limited 
5 image signal) Org into a processed band-limited image signal 
as Proc=N X Org 
[0101] 

In addition, the noise suppressing method described 
above may be carried out by an computer, and a program for 
10 instructing the computer to carry out the noise suppressing 
method may be stored in the computer-readable storage medium 
and such a computer-readable storage medium may be provided. 
[Brief Description of the Drawings] 
[FIG. 1] 

15 FIG . 1 is a schematic block diagram illustrating a construction 
of a noise suppressing apparatus according to one embodiment 
of the invention. 
[FIG. 2] 

FIG. 2 is a flow diagram illustrating a processing sequence of 
2 0 the noise suppressing apparatus. 
[FIG. 3] 

FIG. 3 is a block diagram illustrating an outline of a 
construction of the band-limited-image-signal generation 
unit . 
25 [FIG. 4] 

FIG. 4 is a diagram schematically illustrating the operations 
of generating band-limited image signals. 
[FIG. 5] 

FIG. 5 is a block diagram illustrating details of a general 
30 configuration of the noise suppressing apparatus. 
[FIG. 6] 

FIG. 6 shows conceptual diagrams for illustrating the 
double -angle representation . 
[FIG. 7] 

35 FIG. 7 is a diagram illustrating the definition of the D-A 
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representation . 
[FIG. 8] 

FIG . 8 shows diagrams respectively illustrating four 
two-dimensional filters . 
5 [FIG. 9] 

FIG. 9 is a diagram illustrating a relationship between a 
Laplacian signal and an output of a first-derivative filter 
(the absolute value of the first derivative) . 
[FIG. 10] 

10 FIG. 8 shows views respectively illustrating isotropic 
two-dimensional filters . 
[FIG. 11] 

FIG. 11 shows views respectively illustrating recognition 
models of a line, a point of intersection, an end point, and 
15 noise based on a degree of edge confidence and an index of pixel 
energy . 

[FIG. 12] 

FIG. 12 shows views respectively illustrating examples of 
anisotropic filters having various characteristics. 
20 [FIG. 13] 

FIG. 13 is a diagram schematically illustrating an operation 
for performing the Laplacian reconstruction. 
[FIG. 14] 

FIG. 14 is a detailed block diagram illustrating a construction 
25 of a noise suppressing apparatus according to another 
embodiment of the invention . 

[Explanation of the Reference Numerals] 
100 noise suppressing apparatus 
1 band-limited-image-signal generation unit 
30 2 index-value obtaining unit 

3 noise suppressing unit 

4 image reconstruction unit 
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[FIG. 3] 
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[FIG. 5] 
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[FIG. 9] 



EDGE SIGNAL 



a b 




LAFL*tr/W SPWL 




OUTPUT SIGNAL OF FILTER 



[FIG. 10] 



RELATIONSHIP BETWEEN LAPLACIAN SIGNAL AND OUTPUT OF 

FIRST- DERIVATIVE FILTER (ABSOLUTE VALUE OF FIRST DERIVATIVE) 




Mm. 



(a): FOR TARGET (b): 
RESOLUTION LEVEL 



FOR ONE-LEVEL LOWE* ( c )* TWO-LEVEL LOITER 
RESOLUTION LEVEL RESOLUTION LEVEL 



[FIG. 11] 

< 



(a): LINE 
C= 1 & E= 1 



< 



(b): LINE/NOISE 
C= 1 & E<< 1 



. INTERSECTION POINT 
c ) : OR END POINT 
C<< 1 A E =f 1 



(d): NOISE 

C<< \ & E<< 1 



[FIG. 12] 




0 ' 







135" 



180" 



[FIG. 13] 



fBi 



© 



fB< 



S5 a 






© 


S3' 


► 


S2 




fB= □ DS5 



43 



[FIG. 14] 





[Name of Document] ABSTRACT 
[Abstract] 

[Objective] A noise suppressing apparatus for suppressing 
noise components included in a radiographic image regardless 
5 of an exposure dose is provided which can suppress the noise 
regardless of an exposure dose. 

[Constitution] A noise suppressing apparatus 100 has a 
band-limited-image-signal generation unit 1 for generating a 
plurality of band-limited image signals Bk respectively 

10 representing a plurality of band-limited images belonging to 
a plurality of different frequency bands, based on the input 
image signal Sin; an index-value obtaining unit 2 for obtaining 
an index value indicating a degree of suppression of the noise 
component , based on information indicating an exposure dose 

15 with which the radiographic image has been produced; a noise 
suppressing unit 3 for processing each of the plurality of 
band-limited image signals Bk so as to suppress noise in each 
of the plurality of band-limited images based on the obtained 
index value; and an image reconstruction unit 4 for 

20 reconstructing a processed image signal Sproc which represents 
a noise-suppressed image, from the band-limited image signals 
fBk which undergoes suppression of the noise components. 
[Selected Figure] FIG. 1 
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